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SIGNIFICANCE  AND  EXPLANATION 


A  variety  of  heat-transfer  problems  involving  a  change  of  phase,  such 
as  melting  of  ice,  fusion  of  metals,  solidification  of  alloys,  involves 
basically  two  facts:  a  diffusion  process  (described  for  example  by  the  heat 
equation)  and  a  description  of  the  progressive  change  of  the  region  where  the 
process  takes  place.  Such  phenomena  can  be  formulated  as  free-boundary  pro¬ 
blems  which  are  often  treated  mathematically  using  variational  inequalities. 
The  heat  equation  describes  many  heat  transfer  problems  only  in  the  crudest 
way. 

In  this  paper  we  propose  to  study  such  processes  by  using  a  non-linear, 
more  refined  model  of  heat-conduction.  We  develop  a  mathematical  theory  for 
a  convenient  class  of  variational  inequalities  and  apply  it  to  study  problems 
involving  a  change  of  phase  (for  example  a  Stefan  problem)  where  the  alter¬ 
native  model  mentioned  above  is  employed. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  authors  of  this  report. 


A  PSEUDO-PARABOLIC  VARIATIONAL  INEQUALITY  AND  STEPAN  PROBLEM 
E.  DiBenedetto  and  R.  E.  Showalter^ 

JL  Introduction. 

Let  A  and  B  be  (possibly  multi-valued)  maximal  monotone  operators  and  let 
C  be  a  non-empty  closed  convex  set  in  the  real  Hilbert  space  V.  We  shall  give 
existence  and  uniqueness  results  for  evolution  inequalities  (formally)  of  the  form 

(1.1. a)  u(t)  e  C:  (Au(t))  +  Bu(t)-f  (t),  v-u(t)j  v  >0,  v  e  C  ,  0<t<T, 

(1.1. b)  (Au(0)-v0  ,  v-u(0))v  >  0  ,  v  e  C  , 

2 

where  f  e  L  (0,T;V)  and  v^  e  A(u^)  are  given.  In  Section  2  we  introduce  a  new 

notion  of  weak  solution  of  (1.1)  and  verify  uniqueness  when  A  is  linear  self-cd- 

joint  and  B  is  strictly  monotone.  Existence  of  a  weak  solution  is  proved  in 

Section  3  when  A  is  a  (single-valued)  function  of  the  form  "identity  plus  compact 
* 

operator",  B  is  bounded,  and  A  or  B  is  a  subgradient. 

Variational  inequalities  of  the  form  (1.1)  are  of  interest  on  their  own  as 
extensions  of  corresponding  evolution  equations  of  Sobolev  type  (where  C *  V) .  Early 
work  on  such  inequalities  is  described  in  [2];  we  mention  [6]  specifically  as  a 
source  of  examples  of  initial-boundary  value  problems  for  the  pseudo-parabolic 
partial  differential  equation 

(1.2)  J^(u-aAu)  =  kAu 

with  a  >  0,  k  >  0.  Such  equations  arise  as  models  for  diffusion,  and  they  provide 
an  interesting  alternative  to  the  classical  diffusion  equation  wherein  a»0.  In 
Section  4  we  give  an  example  of  an  initial -boundary-value  problem  consisting  of  a 
highly  nonlinear  partial  differential  equation  of  pseudo-parabolic  type  whose 
solution  is  subject  to  unilateral  constraints.  Existence  and  uniqueness  results  for 
weak  solutions  follow  from  our  abstract  results  on  (1.1). 

^ Department  of  Mathematics,  RLM  8.100,  The  University  of  Texas,  Austin,  TX  78712 
Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 . 


A  one-phase  free-boundary  problem  of  Stefan  type  for  the  equation  (1.2)  Is 
shown  in  Section  5  to  lead  to  the  variational  inequality  (1.1).  This  development 
is  parallel  to  that  of  the  classical  case  a-0  which  is  described,  e.g.,  in  [7]. 
The  existence  of  a  classical  solution  of  a  Stefan  problem  for  (1.2)  in  one  spatial 
dimension  was  given  in  [9]  by  entirely  different  methods. 


2_.  The  Variational  Inequality. 

2 

We  denote  by  L  (0,T;V)  the  Hilbert  space  of  (Bochner)  square-integrable 
functions  on  the  interval  (0,T)  with  values  in  the  Hilbert  6pace  V.  Let 

H*(0,T;V)  denote  the  absolutely  continuous  V-valued  functions  v  whose  deriv- 

dv  2  ^ 

atives  jj:  belong  to  L  (0,T;V).  Denote  the  dual  of  V  by  V  and  recall  the 

2  o  * 

natural  identification  L  (0,T;V)*L  (0,T;V)  ;  thus  we  obtain  the  (dual)  iden¬ 
tification  L2(0,T;V)  c — ►  H*(0,T;V)*  by  restriction.  The  derivative 
HI  9 

*rr  :  Hl(0,T;V)c — »•  L  (0,T;V)  is  a  bounded  linear  operator  which  determines  the  dual 
at 

operator  t  L2 (0* T;V)  H*(0,T;V)  by  the  formula 

T 

(Lf,v>  s  -f  (f(t),v'(t))vdt  ,  f  e  L2(0,T;V)  ,  v  e  H^O.TjV)  . 

°0  v 

The  restriction  of  Lf  to  V-valued  test  functions  is  the  (distribution)  derivative 

Moreover,  for  f  e  H^(0,  T;V)  we  have 
at 

(Lf,v)  -  l~,A  ,  +  (f(0),v(0))  -  (f(T),v(T))  ,  v  e  H^O.TjV)  . 

'<IC  ’u  (0,T;V)  V  V 

Thus,  wc  can  regard  "Lf  +  fff)”  as  formally  equivalent  to  the  Cauchy  operator 

•  I  £  /AN  II 


3. 


We  shall  use  basic  material  on  maximal  monotone  operators  [1] .  Specifically, 
recall  A  cz  VxV  is  monotone  if  [x^  ,  y^]  e  A  for  j*l  and  2  imply 
^xl”x2  '  ^1*^2^V  -  an^  strictly  monotone  if  in  addition  equality  holds  only  if 
x^  •  •  If  (p:  V  >EU  { +“}  is  proper,  convex  and  lower  semicontinuous,  its 

subgradient  defined  by 


B<p(x)  s  {u  e  V:  (u,y-x)y  <  <<>(y)-<p(x)  for  all  y  e  V} 


for  x  e  V  is  maximal  monotone.  More  specifically,  if  C  Is  a  non-empty,  convex 
and  closed  set  in  V,  its  indicator  function 


Ic(x)  h  < 


0 

+  00 


,  x  e  C 
,  x  4  C 


is  prop.er,  convex  and  lower  semicontinuous,  and  we  have  u  e  dl_(x)  if  and  only  if 

V  N  .  V 

x  e  C:  (u,y-x)  <  0  for  all  y  e  C  . 

Thus,  the  subgradient  of  the  indicator  function  provides  a  convenient  method  of 

expressing  the  variational  inequality. 

Suppose  we  are  given  the  pair  A,B  of  maximal  monotone  operators  on  the 

2 

Hilbert  space  V,  a  closed  convex  subset  C  of  V,  f  c  L  (0, T;V)  and  a  pair 

0 

IUq  ,  VqJ  e  A.  Then  a  function  u  is  called  a  strong  solution  of  (1.1)  if  there 
Is  a  pair  of  functions  v,w  such  that 

u, v  e  H^(0, T;V)  ;  w  e  L2(0,T;V)  , 

(2.1. a)  u (t)  e  C:  |^'Uv(t).f(t),  x-u(t)jy  >  0  ,  x  e  C  , 

v(t)  e  A(u(t))  and  w(t)  e  B(u(t))  for  a.e.  t  e  [0,T],  and 

(2.1.b)  (v (0) -vQ  ,  x-u(0))v  >  0  ,  x  e  C  . 


4. 


Note  that  since  u  and  v  are  continuous,  C  is  closed  in  V  and  A  is  closed 
in  V  x  V,  it  follows  that  the  inclusions  u(t)  e  C  and  v(t)  e  A(u(t))  hold  for 
all  t  e  10, T].  Also,  (2.1)  can  be  restated  as 

(2. 2. a)  +  w(t)  +  dlc(u(t))  a  f(t)  , 

v(t)  e  A(u(t))  and  w(t)  e  B(u(t))  for  a.e.  t  e  [ 0, T]  ,  and 
(2.2.b)  v(0)  +  aic(u(0))  a  vQ 

in  terms  of  the  indicator  function. 

We  shall  use  a  weak  notion  of  solution  in  which  it  is  not  required  that 
v  e  H^O.TjV).  Set  K  =  ( u  e  H1  (0,T;V) :  u(t)  e  C,  0  <  t  <  T]  .  Define  a  weak  solution 
of  (1'.  1.)  to  be  a  function  u  for  which  there  is  a  pair  of  functions  v,w  satisfying 

u  e  K  ;  v,w  e  L2(0,T;V)  , 

v (t)  c  A(u(t))  ,  w(t)  e  B(u(t))  ,  a.e.  t  e  [0,TJ  f 

and  for  some  |  c  A(u(T))  we  have 

(2.3)  (Lv+w-f,rj-u>  +  (|,tj(T)-u(T))v  >  (vQ,  t)(0)-u(0))v,  T  e  K  . 

Note  that  if  u  is  a  strong  solution  then  it  is  a  weak  solution  with  |«v(T). 
Moreover  we  have  the  following  elementary  result. 

Theorem  1.  Let  A  be  continuous,  linear,  self-adjoint  and  monotone;  let  B  be 
strictly  monotone.  Then  there  is  at  most  one  weak  solution. 

Proof:  Let  u^  and  u2  be  weak  solutions  and  let  v^  ,  w^^  and  v2  ,  w2  be  the 
corresponding  selections  from  A(u^),  B(u^),  etc.  By  our  assumptions  on  A  we 
have  (after  modification  on  a  null  set)  Vj-A(u^)  e  H1(0,T;V)  and  |j*A(u^(T)) 
for  J“l,2.  Thus  we  have 


<LAuL +v1-f,u2-u1>  +  (Au1(T),u2(T)-u1(T))v  >  (vQ  ,  u2(0)-u1(0))v 
(LAu2 +w2-f,u1-u2)  +  (Au2(T),u1(T)-u2(T))v  >  (vQ  ,  u1(0)-u2(0))v 


For  any  u  £  H  (0,T;V)  we  have 


<LAu,u)  =  l/2((Au(0),u(0))v-(Au(T),u(T))v)  , 
so  adding  the  two  inequalities  and  applying  this  identity  with  u  =  -  u2  gives 


(wrw2  ,  2 


L  (0,  T  ;V) 


+  l/2(Au(T),u(T))v  +  l/2(Au(0),u(0))v  <  0 


Strict  monotonicity  of  B  shows  u^  =  u2  . 

Remarks .  Without  additional  assumptions  we  should  not  expect  uniqueness  of  the 

selectipns  v,w.  For  example,  in  the  extreme  case  C={0),  (2.3)  is  vacuous  and 

2 

we  need  only  choose  v,w  e  L  (0,T;V)  with  v(t)  e  A(0)  and  w(t)  e  B(0)  to  ob¬ 
tain  a  weak  solution.  At  the  other  extreme,  C  =  V,  any  weak  solution  gives  a 

dv  2 

strong  solution  of  the  equation  -^-  +  w=f  in  I*  (0,  T;V)  with  v(0)«=Vq  .  Even 
for  equations,  the  current  uniqueness  proofs  require,  e.g.,  A  or  B  to  be  linear 
self-adjoint.  See  [5]. 


5.*  Existence  of  £  Weak  Solution. 

Our  objective  is  to  prove  the  following  result  on  the  existence  of  weak  solutions 
of  (1.1).  Note  that  each  of  our  hypotheses  concerns  only  one  of  the  three  sources  of 
nonlinearity  in  the  problem;  we  have  not  placed. any  "compatibility*'  conditions  on  the 
operators  A,B  or  the  set  C. 

Theorem  2.  Let  C  be  a  non-empty,  closed  and  convex  subset  of  the  Hilbert  space  V. 
Let  A  and  B  be  maximal  monotone  operators  on  V  and  assume  the  following: 


6. 


(1)  A  is  a  (single-valued)  function  which  maps  bounded  sets  in  V  into  bounded 
sets  in  U,  where  U  is  a  Banach  space  compactly  imbedded  in  V. 

(ii)  B  maps  bounded  sets  in  V  into  bounded  sets  in  V. 

(iii)  Either  A  *  bq>  or  B  =  ckp,  where  <p:  V  *  is  a  convex  and  lower -semi con¬ 
tinuous  function. 

2 

Then  for  each  e  C  and  f  e  L  (0, T;V)  there  is  at  least  one  pair  u,w  such 
that 

(3.1. a)  u  e  K  ,  we  L2(0,T;V)  ,  w(t)  e  B(u(t))  ,  0  <  t  <  T  , 

(3.1.b)  (L(u  +A(u))  +w-f, t)-u)  +  (u(T)  +A(u(T)),rj(T)-u(T))v  >  (uQ +A(u()),-)(0)-u(0))v 
for  q  e  K  , 

and  u  (0)  *  uQ  . 

Remarks  If  in  addition  we  had  A(u)  e  H^(0,T;V),  then  u  would  be  a  strong 
solution  of 

(3.2)  ~(u(t) +A(u(t))) +w(t)+3lc(u(t))  3  f(t)  ,  a.e.  te  10, T] 

This  is  (2. 2. a)  with  A  replaced  by  A  +  I. 

Since  we  do  have  u  e  H*(0,T;V)  and  u(0)  =  Uq,  it  follows  that  (3.1.b)  is 
equivalent  to 

(3.3)  <^  +  LA(u)+w-f,1-u)  +  (A(u(T)),i(T)-u(T))v  >  (A (uQ) , q (0) -u (0) )y  ,  i  e  K  . 

Proof:  We  shall  prove  Theorem  2  in  the  following  steps.  First  we  approximate  (3.2) 
(and,  hence,  (3.1))  by  replacing  by  its  Lipschitz-continuous  Yoshida  approxi¬ 

mation  dl£  ,  c  >  0;  the  resulting  equation  has  a  solution  u  by  [  5  ].  Then  we 

v  C 


7. 


establish  estimates  on  {u  ),  deduce  the  existence  of  a  weak  limit  uslim(u  ), 

£  c 

and  finally  show  u  is  a  weak  solution  of  (3.1). 

The  approx imat ion.  As  an  approximation  of  the  indicator  function  I  we  take 
l£(x)  s  (2e)-1||x-Pc(x)||^  ,  e  >  0  ,  x  e  V  . 

Its  Frechet  derivative  is  oI^(x)  =  e  ^(x-P  (x)),  where  P  is  the  orthogonal  pro- 

v»  c  c 

jection  onto  C,  and  it  is  monotone  and  Lipschitz  continuous.  Thus  B  +  dl^  is 

maximal  monotone  and  we  obtain  from  [  5  ]  the  existence  of  a  pair  u  e  H^(0,T;V), 

2 

W£  e  L  (0,T;V)  for  each  e  >  0  satisfying 


jjr(ue(t)  +  A(u£(t)))  +  w£(t)  +  dl*(u£  (t))  =  f (t)  , 

(3-4) 

we(t)  e  B(u£(t))  ,  a . e.  t  e  [0,T]  , 
and  ue(0)*uq.  This  approximation  (3.4)  is  strongly  suggested  by  (3.2). 

The  estimates.  Consider  the  two  cases  in  (iii) . 

Case  A  =  d<p:  Take  the  scalar  product  of  (3.4)  with  u£  and  integrate;  this  gives 

l/2||ue(t)||J  +  9*(A(ue(t)))  +  r  (vc  ,  ue)v  4-  J*  (dlce(u£),u£)v 

0  0 

■  l/2||u0l|y  +  <p*(A(u0))  +  J  (f,U£)  ,  o  <  t  <  T  , 

•ff 

where  9  (x) =  sup{ (x,y)y-9(y) :  y  e  V)  is  the  convex  conjugate  of  9  [1,  p .  41 ) . 

Since  A  is  bounded,  its  domain  is  all  of  V  so  9(0)  <  «>.  Thus,  we  may  take 
* 

9(0) *  0  and  9  (x)  >  0,  x  e  V,  with  no  loss  of  generality.  Since  B  is  monotone 
(w£(t),u£(t))v  >  (B°(0),u£(t))  for  t  t  [0,T],  where,  e.g.,  B°(0)  e  B(0)  is  the 
minimal  section  of  3  at  0.  Finally  ve  may  assume  0  e  C  arid  thus 


8. 


(Sl*(uc(t)),uc(t))v  >  I*(u£(t))  >  0  ,  0  <  t  <  T  , 

from  the  definition  of  the  subgradient.  These  observations  and  the  preceding 
estimate  give 


(l/*)||u  (t)|£  <  (l/2)||u0||J+ip*(Au0)  +  (||f||  ,  +T||B°(0)|j  )||u  || 

c  V  U  V  U  l/(0,T;V)  V  C  L  (0, T;V) 


0  <  t  <  T  . 


This  implies  that  ||u  ||  is  bounded,  and  from  (i)  and  (ii)  it  follows 

e  L°°(0,  T;V) 

|Jw  ||  and  ||A(u  )  ||  are  bounded  uniformly  in  e  >  0.  Next,  we 

C  L°°(0,T;V)  e  L°°(0,T;U) 

take  the  scalar  product  of  (3.4)  with  u^(t)  and  integrate;  this  gives 


L.^'X  +  I0  (d7  A(ue>  ’  if ) 


V  +  Ic(ue(t» 


<  I?(u  )  +  (||wj  +S|f!|  )|tt'||  2 

L  U  6  LZ(0,T;V)  L  (0,T;V)  e  L  (0,T;V) 


The  monotonicity  of  A  implies  Jthe  second  term  above  is  non-negative  so  we  deduce 

that  ||u'||  o  and  ||lp(u  )||  are  bounded  uniformly  in  e  >  0. 

1*  (0,T;V)  C  e  L“ 

Case  B  =  d<o;  Take  the  scalar  product  of  (3.4)  with  u^(t);  this  gives 


+  ".'4 


(wt(t)+^(ue(t)),u;(t))v 


(f(t),u;(t))v 


The  second  term  is  non-negative  because  A  is  monotone.  The  third  term  is  the 
derivative  of  «p(u  (t) )  +  I«(u  (t))  by  the  chain  rule  [1,  p.73].  Thus  we  integrate 
this  identity  and  obtain 


+  <p(u£(t))  +l£(U(.(t))  <  9(uq)  +  ||f||  2 


L  (0,T;V) 


e"  2 

l/(0,T;V) 


We  can  add  to  B*d<p  a  constant,  by  addtng  the  same  to  f(t),  so  we  may  add  an 
affine  function  to  <p  with  no  loss  of  generality  and  thereby  obtain  «p(x)  >  0 


for  all  x  e  V.  The  preceding  estimate  gives  uniform  bounds  on  |[u'||  9 

e  l/(0,T;V) 

and  ||i*(u  )||  .  Similar  bounds  follow  immediately  on  ||u  ||  and  by  (i) 

C  6  L“  6  L°°  (0,  T  ;V) 

and  (ii)  on  ||A(u  )||  and  ||w  ||  ,  respectively. 

e  L  (0,T;U)  L  (0,T;V) 

The  limit.  From  the  estimates  obtained  above  it  follows  there  is  a  subnet  of  (ue) 
(which  we  denote  again  by  (u^))  for  which 

w-lim(u  )  “  u  in  H*(0,T;V)  ,  and 
£ 

w  -  lim(wc)  =  w  in  L^(0,T;V)  , 
where  "w  -  lim"  denotes  the  weak  limit. 

Lemma  1.  u  e  Ks(v  e  H^O.TjV):  v(t)  e  C,  all  t  e  [0,T]),  u(0)*uQ  and 

w  -  lim(u  (t))*u(t)  in  V  for  every  t  e  [0,T]. 

£ 

Proof:  Let  t  >  0.  For  each  x  e  V  we  have 


(ue(t)-u(t),x)v  =  J  (<ue-u)',x)v  +  (Uq_u  (0), x)^ 


convergent  to  (uq-u(0),x)v  .  By  bounded  convergence 


lim  j  (uc-u,x)y  =  J  (u0-u(0),x)v  =  0 


so  u(0)=un  and  w  -  lim(u  (t) )  =  u(t) .  Next  define  z  (t)  ■  P  (u  (t)),  the 

U  £  C  u  C 

2 

orthogonal  projection  onto  C.  Then  [z^)  is  bounded  in  L  (0,T;V)  so  there 

2 

is  a  subnet  [z. }  which  converges  weakly  to  z  in  L  (0,T;V).  Note  that 
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Iluc(t)-Zc(t)llj  =  2c  I^(ue(t))  <  (const. )c  , 

2 

so  w  -  lim(u^)  =  u  =  z.  Since  the  set  (v  e  L  (0,T;V):  v(.t)  e  C,  a.e.  t  e  C)  contains 
each  and  is  weakly  closed,  it  also  contains  u=z.  Finally  u  e  K  follows  since 

C  is  closed  and  u  is  continuous. 

2 

Lemma  2.  We  have  the  (strong)  limits  lim  A(ue)=A(u)  in  L  (0,T;V)  and 
lira  A(ut(t))  =A(u(t))  in  V  for  every  t  e  [0,T]. 

Proof:  Let  t  e  [0,T].  Since  (A(uE(t)))  is  in  a  compact  set  in  V  there  is  a  sub¬ 
net  (A(u^(t)))  which  converges  (strongly)  to  v(t)  in  V.  But  w  -  lim(u^ (t) ) * u(t) 
in  V  and  A  is  maximal  monotone  so  v(t) =  A(u(t)) .  The  above  applies  as  well  to 

any  subnet  of  (A(ue(t))},  so  the  entire  net  converges  to  A(u(t)).  The  convergence 
2 

in  L^OjTjV)  of  {A(u  ))  to  A(u)  follows  by  the  bounded  convergence  theorem. 

£ 

Lemma  3.  w  e  B(u)  and  lim(w  ,  u  )  7  *  (w,u)  , 

c  LZ(0,T;V)  L  (0,T;V) 

Proof:  It  suffices  to  show  that  [1,  p.  27] 


(3.5) 


lim  sup (w  ,  u  )  -  <  (w,u)  , 

6  LZ(0,T;V)  L  l0, T;V) 


Take  the  scalar  product  of  (3. A)  with  u^  -  u  and  integrate.  From  the  estimate 


Ol*(u  ),u  -u)  >  l!(u  )  -  I  (u)  =  I*(u  )  >  0 

C  e  6  L2(0,T;V)  C  c  C  C  e 


we  obtain 
(w 


t  ■  £  (uc  •  “>L.  +  («(ue  +A(ut>)>"-uc)L2  +  <f>  V”>L* 


11. 


By  taking  the  upper  limit  ve  have 


11m  sup(w£  ,  uf)  2  <  (w,u)  j  +11™  SUP  (^(ue  +  A(ue))*u-ue)12 
L  L  '  /L 


Thus  it  suffices  to  show  the  last  term  is  non-positive.  From  the  identity 


T  1  2  -T 

\\£-  ll  tl^v 


J0  (ue  ’  Vu)v  “  1/20luc(T)l!v“lluOllv)  -J0  K  ’  u)v 


and  Leimna  1  there  follows 


.T  ,  2  -T 

it  ii 


11m  Inf  j  (u^.  ,  ue-u)y  >  l/2(|ju(T)l|J-||u0|p  -  (u',u)v  *=  0  . 


Similarly  from  Lemma  2  and 


-X  x 

J  (A(ue)',ue-u)v=  (Aue(T),ue(T)-u(T))v-J  (A(ue),u;-u')v 


it  follows  that  lim  T  (A(u  )',u  -u)  ■  0. 

J0  c  e  V 


The  solution.  To  show  that  u  is  a  weak  solution.  It  suffices  by  Lemma  1  to  verify 
(3.1.b).  For  any  i)  e  K  it  follows  from  (3.4)  that 


(3.6)  (u'.n-u^  +  (A(u  _)',n-u  )  2  +(w-f,n-u)  2 

e  c  L  <0,T;V)  e  L  (0,T;V)  C  C  L  (0,T;V) 


-  (di:(u  ),u  -n) 

C  c  c  h  (0,T;V) 


From  the  definition  of  subgradient  it  follows  that  the  right  side  is  greater  than 


je(u  )  -  l*0i) *  I*(u  ),  hence,  non-negative.  Consider  the  first  term  on  the  left 
C  c  C  C  c 


side.  We  have  from  w  -  lim  u  (T)*u(T)  that 
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<u\u)  2-  (1/2)(!!u(T)||^-|!u0I!2)  <  lint  inf  (1/2)  (||ue(T) ||J-||tt0llJ) 
L 

*  lim  inf  (u^.  ,  u£)  2  , 


so  there  follows 


<Lu,  tj-u)  +  (u(T),r)(T)  -u(T)). 


(uQ  ,  r,(0)-u0)v=  (u',tj-u)  2  >  11®  sup(u^  ,  n-u^  2 

L  L 


Concerning  the  second  term,  we  obtain  from  Lemma  2 


<L(Au),n-u>  +  (A(u(T)),n(T)-u(T))v  -  (A(u0),n(0)-u0)v 

-  lim(<L(Aue),t,-ue)  +  (A(oe(T)),T)Cr)-ue(T))v-(A(u0),r,(0)-u0)v) 

-  lim(A(u  )',n-u  )  , 

6  LZ(0,T;V) 

Finally,  Lenina  3  identifies  the  limit  of  the  third  term,  so  by  taking  the  "lim  sup" 
in  (3.6)  we  obtain  (3.1.b). 


( i..  A  Pseudo-Parabolic  Inequality. 

When  our  results  from  above  are  used  to  describe  initial -bound ary -value  problems 
for  partial  differential  equations  or  inequalities,  it  is  usually  more  convenient  to 
express  them  in  terms  of  the  equivalent  notion  of  a  maximal  monotone  operator  Q 
from  the  Hilbert  V  to  its  dual  V*.  Thus,  letting  9t:  V  V  be  the  Riesz  iso¬ 
morphism  given  by  the  scalar  product, 

S*x(y)  -  (x,y)y  ,  x,y  e  V  , 

we  say  flcVxV*  is  monotone  if  the  composite  operator  As 91  ^  •  fl  is  monotone  in 

V  x  V  and  maximal  monotone  if,  in  addition,  Rg(£R  +  fl)»V*.  We  can  easily  state 

Theorem  2,  for  example,  in  this  context.  Thus,  we  are  given  a  set  C  closed  and 

★ 

convex  in  V;  3  and  B  are  maximal  monotone  operators  from  V  to  V  satisfying 
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hypotheses  corresponding  to  (i),  (ii)  and  (iii).  Then  for  each  e  C  and 

2  iHr 

f  e  L  (0,T;V  )  there  is  a  pair  of  functions  u,w  satisfying 

(4.1. a)  u  e  K  =  (v  e  H^O.TjV):  v(t)  e  C  for  0  <  t  <  T)  , 

(4.1.b)  w  e  L2(0,T;V*)  ,  w(t)  e  S(u(t))  for  a.e.  t  e  [0,T]  , 

(4.1.c)  «(u)  e  L2(0,T;V*)  and  u(0)  =  , 

(4.1.d)  <L(£R  +  fl)(u) +w-f,n-u>  +  (£  +  G)u(T)0l(T)-u(T)) 

>  (3M-fl)uQ(ri  (0)  -u(0))  ,  for  T)  e  K  . 

2*1  * 

In  this  setting  the  linear  operator  L:  L  (0,T;V  )  H  (0,  T;V)  is  given  by 

T 

<Lg,u>  -  -f  g(t)(u'(t))dt  ,  g  e  L2(0,T;V*)  ,  u  e  H^O.TjV)  . 
u 

Since  u  e  H*(0,T;V),  the  inequality  (4.1.d)  is  equivalent  to 

<~3hi+LG(u)  +w-f,r)-u> +fl(u(T))(q(T)-u(T))  >  fl (uQ)  (q (0) -u (0) )  , 
for  i)  t  K  . 

We  shall  describe  an  example  of  a  partial  differential  equation  of  pseudo¬ 
parabolic  type  which  is  to  be  resolved  subject  to  unilateral  boundary  constraints. 
(A  similar  equation  with  constraint  over  the  entire  region  will  be  given  in  the 
next  section.)  Let  G  be  a  bounded  open  set  in  Rn  which  lies  on  one  side  of 
its  boundary  dG;  assume  dG  consists  of  two  disjoint  parts  Tq  and  r,  and 

let  n(s)  *  (ti^ (s), . .  .,nn(s))  be  the  unit  outward  normal  at  each  point  s  c  dG. 

1  2 
H  (G)  is  the  Sobolev  space  of  those  v  e  L  (G)  for  which  all  derivatives 
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c3v  2  1 

*  DjV,  1  <  J  <  n,  belong  to  L  (G);  we  set  DQv  =  v.  Let  V  =  {v  e  H  (G): 

v |  ™  0} ;  by  v|  we  mean  the  trace  of  v  on  rft  (see  [8,10]).  For  v  c  V 

*0  l0  u 

2  00 
we  denote  by  y(v)  e  L  (r)  the  trace  of  v  on  r.  Let  r^  e  L°° (G)  and 

00 

r^  e  L  (r)  be  non-negative  and  define 


(4.2)  £Ru(v)  =  (u,v)v  =  J  |  I  DjUDjV +  rQuvj  +  J  r.jY(u)Y(v),  u. 


v  e  V 


It  follows  by  a  compactness  argument  that  (4.2)  is  equivalent  to  the  usual  li  (G) 
scalar-product  if  any  one  of  Tq  or  (x:  r^(x)  >  0)  or  [s:  r^(s)  >  0}  has 
strictly  positive  measure,  and  we  assume  this  hereafter. 

The  operator  <2  is  given  by  a  pair  of  continuous  (maximal)  monotone  functions 
O.Q  ,  a^:  R  ->-R  which  are  linearly  bounded: 


|a.j(z)  |  <  Q(1  +  |*|)  ,  zeR,  j  «  0, -1 


for  some  constant  Q  >  0.  We  define 


(4.3) 


fl(u)(v)  *  J  oQ(u)v  +  J  a_1(Yu)yv  ,  u,v  e  V 


This  operator  is  a  subgradient  (in  fact,  a  Gateux  derivative)  and  is  bounded  from 
V  to  UsL^(G)  xL^(r).  Since  the  imbeddings  V  c— ►  U  and  U  =  U*  c-*  V*  are 
compact,  the  hypothesis  (i)  is  fulfilled. 

The  operator  B  will  be  specified  by  a  family  of  maximal  monotone  operators 
(J^:  R  *  R  which  are  linearly  bounded: 

j w |  <  Q(1  +  |z|)  for  w  e  Pj(z)  ,  z  e  R  ,  -1  <  k  <  n 
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for  some  Q  >  0.  Each  P°r  a  corresponding  convex  continuous  ^ 

We  then  define  S  =  Sep  where 


R  *  R. 


o  n  f* 

<p(v)  a  \  £  94(D4v)  +  V 

JG  j-0  j  j  Jr 


«P_1(yv) 


v  e  V  . 


Thus  $  is  given  (formally)  by 

(4.4)  »(u)  «  £  D*(3.(D.u)  +y*P.1(Yu)  ,  u  e  V  . 

k*0 

To  be  precise,  we  have  F  e  B(u)  if  and  only  if  there  exist  f^  e  p^(D^u)  in 
L2(G),  0  <  1:  ,  and  f  .  e  P  ,  (yu)  in  L2(r)  for  which 

—I  -i 

n  o 

F(v)  -  J  £  fk\v  +  J  f _! (Yv)  ,  v  e  V  . 

G  k®0  P 

By  restricting  each  of  the  functionals  SR(u),  fl(u)  and  33 (u)  to  test  functions 

00 

Cq(G)  we  obtain  the  corresponding  distributions  over  G 


(4. 5. a) 
(4.5.b) 
(4.5.c) 


*0(«) 

V“> 

Vu) 


-&»  +  rou  • 

a0(u)  , 
n 

-ED.p.CD.u)  +  p0(u) 
j«l  3  3  3  u 


u  e  V  , 


where  the  multi-valued  $q  is  interpreted  as  before.  The  respective  differences 
are  given  (by  Green's  theorem  for  the  first  and  last  cases)  for  sufficiently  regular 
u  by 


j 
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(4. 6. a)  9tu(v)  -  \  9L(u)v  ■ 

JG  0 

(4.6.b)  fi(u)(v)  - J  aQ(\i)v  -  J  (yu)yv  , 

(4.6.c)  S(u)(v)  -  SQ(u)v  ■  J  ^  I  pj(DjU)nj  +  P_1(yu)J  yv  ,  v  c  V  . 

Thus  we  have  realised  the  operators  a,  fl  and  8  as  the  sum  of  a  distribution 

over  G  (4.5)  and  a  boundary  part  over  T  (4.6).  See  [5]  for  details. 

The  remaining  data  is  given  as  follows.  Let  C«{v  e  V:  y(v)  >  0  a.e.  on  r) 

and  let  uQ  e  C  be  specified.  Suppose  FQ  e  L2(Gx  [0,T])  and  gQ  e  L2(rx[0,T]) 

2  ★ 

are  given  and  define  f  c  L  (0,T;V  )  by 


f(t)(v)«j*  FQ(.,t)v+J  g0(-,t)Yv  , 

g  r 


v  e  V  . 


With  the  preceding  data  as  given,  the  solution  u,w  of  (4.1)  is  a  generalized 
solution  of  the  pseudo-parabolic  problem 


(A.  7)  ) 


^<V  +  fl0 <u)  >  +  ®0  <u)  "  F0  ln  G  x  <°’  T> 

u(x,0)  -  uQ(x)  ,  x  e  G  , 
u  ■  0  on  rQx  [0,T] 

u  >  0  ,  A(u)  >  0  ,  A(u)(u)  ■  0  on  Tx  [0,T] 


idiere  A  is  the  boundary  operator  obtained  from  (4.6)  as 


(4.8) 


A(u) 


^(lS  +  r1Y(u)+a-1(Y«)) 


n 

+  S  3,(D  u)n.  +  p  (yu)  . 
J-l  3  3  3  ~l 


The  other  operators  in  (4.7)  are  given  by  (4.5).  Note  specifically  that  the  multi¬ 
valued  operators  are  to  be  interpreted  precisely  as  was  done  above  following  (4.4). 

5.  A  Stefan  Problem. 

We  consider  a  problem  of  heat  diffusion  involving  a  solid-liquid  phase  change 
at  a  prescribed  temperature.  One  application  we  have  in  mind  is  the  melting  of  ice 
(initially  at  temperature  Kero)  suspended  in  a  reservoir  or  porous  medium.  The 
novelty  in  this  treatment  is  that  we  assume  the  heat  diffusion  is  governed  by  the 
pair  of  equations 

—  m  kAtp  ,  0  *  q>  -  aAcp  . 

Chen  and  Gurtin  [3]  introduced  such  a  model  for  heat  conduction  in  non-simple 
materials  where  the  energy,  entropy,  heat  flux  and  thermodynamic  temperature  9(x, t) 
depend  on  the  conductive  temperature  <p(x,t)  and  its  first  two  spatial  gradients. 
Here  the  heat  flux  is  determined  by  the  conductive  temperature  and  the  phase  is  de¬ 
termined  by  the  thermodynamic  temperature.  Thus  9  >  0  in  the  region  occupied  by 
water  and  9  =  0  corresponds  to  the  frozen  region. 

We  describe  the  geometry  of  the  problem.  Let  the  bounded  domain  G  in  Rn 
be  the  medium  in  which  the  ice/water  is  suspended  and  let  its  boundary  dG  consist 
of  two  disjoint  pieces,  Tq  and  .  Set  fi  =  Gx  (0,T),  where  T  >  0,  and  note 
that  its  lateral  boundary  is  B^  U  ,  where  B^  =  rjX(0,T)  for  j  =  0, 1.  The 
water-region  *  ( (x, t)  e  0:  9(x,t)  >  0)  is  separated  from  the  ice-region 
r<Q»{(x,t)  e  0:  9(x,t)*0)  by  an  interface  S  which  is  the  phase  boundary.  The 

j.  ^  n 

unit  outward  normal  on  oQ^  is  denoted  by  N«  (N^  ,  Nfc),  Nx  c  R  .  If  V(t)  is 
the  velocity  in  Rn  of  the  interface  at  time  t,  then  it  follows  by  the  chain 
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rule  that  V(t) -N^  ■  0  on  S.  Set  n  =  Nx/ liNx U *  t*'e  unit  outward  normal  in  Rn 

of  the  lateral  boundary  of  .  Of  course  n«Nx  on  ,  and  » 0  where  t«0 
or  t-T. 

The  problem  is  formulated  as  follows.  We  are  given  the  conductivity  k  >  0, 
temperature  discrepancy  a  >  0,  and  latent  heat  b  >  0,  of  the  material  and  a 
constant  h  >  0  representing  conductivity  across  the  lateral  boundary  .  The 
initial  thermodynamic  temperature  9q(x),  x  e  G,  and  applied  conductive  temperature 
g(x, t),  (x,  t)  e  B^  ,  are  given  with  9Q=0  on  ro  ’  °0  >  0  °n  ri  *  an<1  8  >  0. 

The  local  form  of  the  problem  is  to  find  a  pair  of  non-negative  functions  0,tp  on 
G  for  which  we  have 


(5.1. a) 

=  kAcp  ,  9  *  tp  -  aAip 

in 

°1 

<5.1. by 

k  ^  +  bV (t)  *n  **  0 

on 

S 

(5.1.c) 

k  |^  +  h(«p-g)  **  0 

on 

F1 

(5.1.d) 

9  ■  0 

on 

ro 

(S.l.e) 

®(*»  0)  -  9q 

on 

G  . 

Note  that  if  0,<p  is  a  solution  of  (5.1)  and  9q  >  0,  then 
(5.2)  (a/k)  +  0  -  <p  in  n 

so  it  follows  that  <p"0  on  U  S.  Since  g.>  0,  the  maximum  principle  for  the 

elliptic  equation  in  (5.1. a)  on  the  region  G(t)  s  (x  c  G:  (x,t)  e  G^  shows  that 
<p  >  0  in  and  0  on  S.  Thus  <  0  on  S  and  G(t)  is  increasing 


with  t . 


We  shall  show  that  the  problem  (5.1)  leads  to  a  variational  inequality  of  the 

form  (1.1).  Define  V  =  (v  e  H1  (G) *.  v|  =0]  as  before.  Regarding  regularity  of 

0  1 

a  solution,  we  assume  0q  e  V,  0:  (0,T]  V  is  absolutely  continuous,  e  L  (0,T;V), 
and  (c.f.  (5.2)) 

(5.3)  k  ^dt^*+  »  a,e*  t  6  t0,Tl  ’ 

* 

Define  the  continuous  linear  JR:  V  >  V  by 


»u(v)  »  J  k(^i'Vv)dx  +  J  h(uv)ds  ,  u,v  e  V 

c  rx 


For  a  test  function  v  c  Cq((0,T),V)  we  obtain 


T 

r  3Up(t)(v(t))dt  -  J  kV«p-Vv  dxdt  +  J  h«pv  dsdt 


r  (-kA<p)v  dxdt  +  f  k  Vip'H  v  dsdt  +  f  h<pv  dsdt 

\  x  \ 

-Ja  If v + JBihgv + Js  bV 


from  (5.1).  Furthermore  we  have 


Tv -I  fc-J 


'ft,  ft 


in  the  sense  of  V*-valued  distributions,  where  H(s)  *  1  for  s  >  0  and  H(s)*0 
for  s  <  0  is  the  Heaviside  function.  We  can  summarise  the  above  calculations  as 


(5.4) 


(9  +bH(9))  +  Sfcp  =  (gh)j 


In 


L1(0,  T;V*) 


where  we  define 


(hg). 


(t)(v)  *  J  hg(s,  t)v(s)ds  ,  v  €  V  ,  t  e  [0,T] 


Combining  (5.3)  and  (5.4)  we  find  that  the  absolutely  continuous  function 
9:  [0,T]  >  V  satisfies 


(5. 5. a)  £  (9+(a/k)9l(9)+bH(9))+a(9)*  (gh)^  in  ^(O.TjV*), 

(5.5.b)  9(0)  =  9Q  ,  and 


(5.5.c)  9(x,t)  >  0  ,  a.e.  x  e  G  ,  t  e  [0,T]  . 

If  we  integrate  (5. 5. a)  and  follow  the  suggestion  in  [7]  to  set 


0(s)ds 


f(t)  s  (1  +  (a/k)$R  +  bH)90  -b 


(hg)_  (sXds 

ri  \ 


» 


there  follows 


~  (1+  (a/k)30u  +  JRu  -  f(t)  -  b(l-H(9))  . 

Finally  we  note  that  H(u)=H(9)  since  G(t)  is  increasing  in  t,  hence, 
u(l-H(9) )  =  0  in  n. 

The  preceding  computations  show  that  u  e  H*(0,T;V)  and  it  satisfies 
u(0)  -  0, 
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(5. 6. a)  u(t)  >0  in  V 

<5.6.b)  (9t+ (k/a)I)u(t)  +  (k/a)3tu(t)  >  f(t)  in  V*,  and 

(5.6.c)  ^~(£R+(k/a)I)u(t)  +  (k/a)$u(t)-f  (t)J  (u(t))  =  0  ,  0<t<T. 

Setting  C  =  {v  e  V:  v  >  0  a.e.  in  G}  we  see  that  u  is  a  strong  solution  of 

(4.1)  with  G  =  (k/a) I,  S=(k/a)9l,  and  uQ*0.  Theorem  1  asserts  the  uniqueness 

of  a  solution  of  (5.1)  under  conditions  considerably  weaker  than  those  leading  to 

(5.6).  Theorem  2  establishes  the  existence  of  a  weak  solution  with  certain 

1  * 

additional  regularity  properties.  In  particular  flu  e  H  (0,T;V  )  since  fl  is 
continuous  and  linear,  so  (4.1.d)  is  equivalent  to 

(“■(£Ru  +  flu)  +w-f,  q-u)  >  0  ,  for  q  e  K  . 
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